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Abstract 

We consider an electrical chain of N nonlinear capacitors coupled by 
linear inductors assuming that voltage dependence of capacitors represents 
an even function. We prove that only 5 symmetry determined nonlinear 
normal modes (NNM) can exist in the considered system. The stability of 
all these dynamical regimes for different N is studied with the aid of the 
group-theoretical method [Physical Review E 73 (2006) 36216] which allows 
to simplify radically the variational systems appearing in the Floquet stability 
analysis. The scailing of the voltage stability threshold in the thermodynamic 
limit N — y oo is determined for each NNM. 

Keywords: nonlinear dynamics, lattice models, nonlinear normal modes, 
invariant manifolds, group-theoretical methods 

1. Introduction 

In recent years, studying of nonlinear vibrations in mesoscopic systems 
of different physical nature received much attention. In particular, discrete 
breathers and various soliton-like excitations were actively studied in can- 
tilever arrays [H, (2), granular crystals |3|, Josephson junction lattices 0, 0|, 



photonic crystals [6j], electrical lines prill], etc. It is essential that in con- 
trast to crystals, where only indirect experiments are possible, in mesoscopic 
systems one can often observe different dynamical objects directly. 



In Ref. 10], the existence and linear stability of the so-called 7r-mode 



(zone boundary mode) in one-dimensional nonlinear electrical lattice were 
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studied. This system, constructed on a silicon substrate by CMOS technol- 
ogy, represents a chain of nonlinear (voltage-dependent) capacitors coupled 
by linear inductors. This chain is depicted schematically in FigJIJi. Using 
periodic boundary conditions, one can imagine that considered circuit repre- 
sents a ring of N cells. For the case N = 4, such ring is depicted in FigJTb. 

The dependence C(V) of the capacity C on the voltage V can be rather 
complex, but the authors of the discussed reference used the following simple 
form of the function C(V): 

C(V) = C (l-bV 2 ), (1) 

where Co and b are positive constants. 

Applying the Kirchhoff laws for quasi-stationary current to the circuit in 



FigJU one can obtain the following equations (see Ref. (l0|): 



L *tj_y y dQ{V 3 )_ 

Here Q(Vj) = Q{Vj(t)) is a charge on j-th capacitor at instant t [the voltage 
on this capacitor is equal to Vj(t)}, while Ij = Ij(t) represents current through 
the inductor Lj 0. 

Taking into account relations Ij = Qj and Qj(t) = L J C(V)dV = 
CoVj(l — fVj 2 ), one can reduce the two first-order Eqs.fl2]) to the system 
of differential equations of the second order with respect to voltages Vj(t): 

Vj{\ - bVf) - 2bV j (V j ) 2 = 7(^-1 - 2Vj + V j+ i). (3) 

Hereafter all inductors are supposed to be identical: Lj = Lq (j = 1..N) and 
7 = \/{LqCq). The appropriate scaling of variables in Eq.flH]) allows one to 
suppose 6 = 1. Finally, taking into account periodic boundary conditions, we 
obtain the following dynamical model which is studied in the present paper: 



Vj(l - Vf) - 2V J (V j ) 2 = 7(^-1 - 2V5 + V j+1 ), (4a) 
V (t) = V N (t), V N+1 (t) = Vx(t). (4b) 



1 Note that there is no resistance and, therefore, no energy dissipation in the considered 
model ((21). More realistic models can be found in [Tlj ] . 
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The so-called zone boundary mode, or 7r-mode, which was investigated in 
Ref. [10|, can be written in the form 

<t> x = {V(t), -V(t)\V(t), -V(t)\V(t), -V(t)\...}. (5) 

This means that voltages on every pair of neighboring capacitors are opposite 
in sign. Obviously, such dynamical regime can exist only in chains with an 
even number of cells. 

There are many papers devoted to study 7r-mode in monoatomic chains 



of the Fermi-Pasta-Ulam (FPU) type [12l4l8j|. However, the model (jl]) es- 
sentially differs from the standard FPU-equations. In Appendix of Ref. 10], 
the authors using some approximation reduce equations of the considered 
electric chain to those of the FPU-beta model, but we prefer to study the 
problem directly in the form (jl]). 

Stability of 7r-mode depends on the number of chain cells (N). Since this 
mode represents periodic regime, its stability can be investigated with the 
aid of the standard Floquet method. However, dimension of the variational 
system (the system of the equations linearized in the vicinity of the considered 
dynamical regime) and corresponding monodromic matrix is equal to 2N. 
Because of this reason, studying of the linear stability of 7r-mode can be very 
difficult in the case N 3> 1, especially when N — > oo. A special method for 



studying this problem was developed in [10]. However, this method is based 
essentially on the specific structure of the dynamical model of the electrical 
chain, and seems to be rather difficult. 

On the other hand, a general group-theoretical method for splitting (de- 
composing) variational system into subsystems, whose dimensions can be 
considerably smaller than that of the original variational system, was devel- 



oped for dynamical models with discrete symmetry in our papers [181 . 119 
This method is based only on symmetry-related arguments and uses the appa- 
ratus of irreducible representations of the symmetry group of the considered 
dynamical regime. 



In Refs. [la . 1 1 61 ] , we used this method for analyzing linear stability of 



all possible symmetry-determined Rosenberg nonlinear normal modes in the 



FPU-a and FPU-/3 chains. In Refs. [20j, |21|, it was applied for studying 
stability of discrete breathers and quasibreathers in 2D scalar dynamical 
models on the plane square lattices. 

In the present paper, we use the above group-theoretical method for an- 
alyzing stability of the symmetry-determined nonlinear normal modes in the 
electrical model (HI). 
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The paper is organized as follows. The possible symmetry-determined 
nonlinear normal modes (NNMs) in the model (jl]) are discussed in Sec. 2. 
In Sec. 3, we briefly review the group-theoretical method for studying their 
stability. In Sec. 4, we apply this method for splitting variational systems for 
all NNMs permissible in the model (jl]). Sec. 5 is devoted to stability analysis 
of 7r-mode, while stability of others NNMs are discussed in Sees. 6-9. In 
Conclusion, we summarize the results on stability of NNMs in the considered 
electrical chain and discuss perspectives of applying our group-theoretical 
method for studying stability of quasiperiodic nonlinear vibrations (bushes 
of NNMs). 



2. Symmetry-determined nonlinear normal modes 



in 



The concept of nonlinear normal modes was developed by Rosenberg 
22| (more information on these dynamical objects can be found in 23]). 



In the dynamical regime corresponding to a given NNM, vibrations of all 
dynamical variables are described by one and the same time-dependent func- 
tion f(t). For the electrical model (JH), this means that all voltages Vj(t) 
satisfies the relation 

V j (t)=c j f{t),j = l..N, (6) 



where Cj are constant coefficients. 

Note that convenient linear normal modes also satisfy the definition ([6]) 
with f(t) = sin (cut + <p ), where u and ip represent frequency and initial 
phase, respectively. 

Let us emphasize that Rosenberg nonlinear normal modes can exist only 
in some very specific classes of systems, for example, in systems whose poten- 
tial energy is a homogeneous function of all its arguments. However, existence 
of NNMs in systems with general type of interactions can be ensured by the 



presence of some group of discrete symmetry [15l. Il6l. |24H26[|. Hereafter, we 
call such modes symmetry-determined nonlinear normal modes, or simply 
NNMs because only modes of this type are considered in the present paper. 

It is very essential, that there exists only a small number of symmetry- 
determined^ NNMs in any lattice model. We discussed the cause of such situ- 
ation in 



16, 26 



The matter is that the attempt of construction NNM with 
large cell in vibrational state leads to appearance of quasiperiodic dynamical 
regimes which represent bushes of NNMs 25, 27, 28[ (see also 26, 29, 30]). 



Let us give a brief review on this subject. 
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All possible dynamical regimes in physical system described by discrete 
symmetry group Go (the dynamical equations of the system are invariant 
according to all transformations of Go) can be classified by subgroups of the 
group G . 

It is easy to check that equations (jlj) are invariant under the action of 
the group Go which is isomorphic to the symmetry group D n of the FPU-/3 
chain. This group can be described by generators which we define by their 
action on the vector 

V = {V 1 (t),V 2 (t),...,V N (t)} 

determining a given vibrational state of the electrical chain. 

Transitional transformation a induces the cycle transposition of variables 
Vj (note that periodic boundary conditions are assumed): 

aV = {V 2 (t), V 3 (t), ...V N (t), Vrit)}. (7) 

Transformation i generates inversion of the components of the vector V with 
respect to an arbitrary chosen "middle" of the chain: 



W = {V N (t), V N ^(t), V 2 (t), V 1 (t)}, 

Finally, the transformation u changes signs of all variables Vj(t) without their 
transposition: 

uV={-V 1 (t),-V 2 (t),...,-V N (t)}. 

It is easy to check that the action of all generators of the group Go, i.e., 
a, i, u, transforms the dynamical system into itself and, therefore, they 
actually are invariant transformations for the considered electrical chain. 

The orders of the generators a, i, u are equal to N, 2, 2 respectively, i.e. 

a N = e, i 2 = e, u 2 = e, (8) 

where e is the identical element of the group Go- Thus, Go consist of 4iV 
symmetry elements representing all products of the above three generators. 
Note that Go is a nonabelian group (for example, la ^ ai). 

Every subgroup Gj (j ^ 0) of the group G determines a specific dy- 
namical regime that can take place in the model (j4j). The subgroup Gj we 
call the symmetry group of the corresponding regime since the vector V(t), 
describing this regime, is invariant under the action of all elements of Gj. 
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Let us consider the subgroup Gj C Go that consists of identical element 
e and the transformation a 2 for the electrical chain with N = 6. We look 
for the dynamical regime V(t) which is invariant under the action of this 
subgroup: 

a 2 V = V. (9) 

Taking into account that 

a 2 V={V z ,Vi,V h ,V 6 ,Vi,V2), (10) 

we obtain the following relations: 

v 3 = V u V 4 = V 2 , V 5 = V 3 , V 6 = V A , V x = V s , V 2 = V 6 . (11) 

Therefore, the six-dimensional vector V is determined by only two indepen- 
dent variables Vi(t) and V 2 (t): 

V = {V 1 ,V 2 \V 1 ,V 2 \V 1 ,V 2 }. (12) 

With the aid of relations ffTTl) between dynamical variables Vj(t)(j = 1..6) 
we reduce dynamical equations (j4]) to two independent equations with respect 
to variables V\{t) and V 2 (t) (all other equations (j4j) turn out to be equiva- 
lent to them). Thus, we obtain a two-parametric quasiperiodic dynamical 
regime fll2p : there are two basic frequencies and their integer combinations 
in the Fourier spectrum of such type of vibrations. 

However, only periodic vibrations are considered in the present paper and, 
therefore, we must select from all symmetry-determined dynamical regimes 
only those which are one-parametric (only one basic frequency presents in 
their Fourier spectrum). Obviously, these regimes represent Rosenberg non- 
linear normal modes. 

The problem of finding all one-parametric symmetry-determined dynam- 
ical regimes in a monoatomic chain can be solved with the aid of the group- 
theoretic method that was described in details in [l6[ . (Note that the similar 
problem for all periodic structures described by any of 230 space groups was 



solved in 26 



Below, we present only final results of the selection of all symmetry- 
determined NNMs that are possible in chains (jlj) for different number iV of 
electrical cells. 

It was found that only five types of such NNMs can exist in the consid- 
ered electrical chains. Let us emphasize once more that this result does not 
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depends on the specific structure of Eqs. fll]) — it is a consequence only of 
symmetry group of the considered lattice model. 

Below, we list all symmetry-determined NNM for the model fl4]). 



Gi = 


[a 2 , iu] : = 


K(t){l,-l|l,-l|l, 


— 1 ... 1, — l}(7r-mode). 


(13) 


G 2 = 


[a 4 ,iu] : 2 = 


K(t){l,l,-l,-l|l, 


1,-1, -1|...}. 


(14) 


G,= 


[a 4 , ai] : </> 3 = 


V(f){l,0,-l,0|l,0 : 


,-l,0|...}. 


(15) 


G4 = 


[a 3 , iu] : 04 = 


K(t){l,0,-l|l,0, - 




(16) 


G 5 = 


[a 3 u, aiu] : 5 


= V(t){l, 1,0,-1,- 


-1,0|1,1,0,-1,-1,0|...}. 


(17) 



Let us comment on the notations used in Eqs. f[T3"j) - f[T7]) . Each subgroup 
Cj C Gq{j = 1..5) is determined by two generators which are shown in square 
brackets. (Let us remind that the original group Go is determined by three 
generators: Go = [a,i,u]). For example, the generators a 2 and iu correspond 
to the subgroup G\ = [a 2 ,iu]. The action of a 2 on the vector V(£) for the 
chain with N = 6 was shown in Eq. fllOl) . while the second generator, iu acts 
on this vector as follows: 

iuV = {-V 6 (t), -V 5 (t), -Vt{t), -V 3 (t), -V 2 (t), -Vi(t)}. (18) 

Demanding invariance of the vector V(t) under the action of a 2 , we have 
obtained dynamical regime (|12p . The presence of the second generator (iu) in 
the subgroup G\ leads to an additional relation between variables Vi(t), ^(i): 

V 2 (t) = -Vx(t). (19) 

Therefore, we obtain the one-dimensional regime representing 7r-mode ffl3|) . 
Note that in Eqs. (IT5|) - (IT7|) we separate by vertical lines the different cells of 
the vibrational state of the electrical chain. In both cases, (|12|) and (fT3l . this 
cell is twice larger than that of equilibrium state and includes two capacitors. 

Vibrational cells including four neighbor capacitors correspond to the 
modes fl!4p and (fl5|) . and those embrasing three and six capacitors are asso- 
ciated with the modes ([TBI) and (TPTj) . respectively. 

It is very essential that the modes with larger vibrational cells cannot 
exist in the considered electrical chain. Indeed, if we try to construct a 
dynamical regime of such type, it turns out to be quasiperiodic because of 
dependence of more than one independent variable (it represents a bush of 
nonlinear normal modes with dimension m > 1). 
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Substituting Eqs. (IT5|) - (IT7|) into dynamical equations (@J, we obtain only 
one independent equation with respect to the variable V(t): 



V{t)[l - V 2 (t)] - 2V(t)[V(t)] 2 + ftV(t) = 0, (20) 

where the following values of the parameter ft correspond to the NNMs shown 
in Eqs. (IH-dlTD: 

G\ = [a 2 , iu] : ft — 4, 

G 2 = [a 4 , iu] : ft — 2, 

G 3 = [a 4 , at] : ft = 2, 

G4 = [a 3 , iu] : ft = 3, 

G5 = [a 3 u, aiu] : ft = 1. 

Hereafter, we call Eq. (|20|) governing equation for the corresponding non- 
linear normal mode. 



3. The method for studying nonlinear normal mode stability 

A given NNM represents periodic dynamical regime and, therefore, the 
standard Floquet method can be used for studying its stability. Using this 
method, we linearize the original system of nonlinear differential equations 
near this NNM and, as a result, we obtain the so called variational system. 
This is a system of N second-order differential equations with time-periodic 
coefficients, whose period T is equal to that of the given NNM. Then we 
reduce this system to the system F of 2N first-order differential equations 
and construct for it the monodromic matrix. The successive columns of this 
matrix can be obtained by integrating 2N times the system F over period T 
using the columns of 2N x 2N identical matrix as the corresponding initial 
conditions. 

As was already noted, the above-mentioned Floquet procedure is very 
difficult for the case N ^> 1, especially when iV — > 00. For the dynamical 
systems with discrete symmetries these difficulties can be often overcome 



with the aid of the group-theoretical method developed in our papers [18|, |19 
Below, we give brief outline of this method. 

Let us consider the symmetry group G of the dynamical model (0J. 
Every element g G Go induces a certain transposition of the variables Vj (j = 
1..N), possibly, accompanied by changing signs of some of these variables. 
Noting that Vj and Vj are transformed under the action of the element g as 
variables Vj, we conclude that this element leads to the same transposition of 
equations (jl]) with the same changing of their signs as the set Vj (j = 1..N). 
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Obviously, arbitrary chosen transposition of variables Vj may not belong 
to the group G . For example, if we transpose only two variables, V 2 <-> V 3 , 
without any transposition of other variables, equations (Jl]) change essentially. 
Indeed, equations 2 and 3 of the system (TJ]) read 

V 2 (l - V 2 ) - 2V 2 (V 2 ) 2 = 7G4 - 2V 2 + V 3 ), 
V 3 (l - V 3 2 ) - 2V 3 (V 3 ) 2 = 7 (V 2 - 2V 3 + V 4 ). 

After the above transposition these equations acquire the following form: 

V 3 (l - V 3 2 ) - 2V 3 (V 3 ) 2 = 7(^1 - 2V 3 + V 2 ), 
V 2 (l - V 2 ) - 2V 2 (V 3 ) 2 = 7(^3 - 2V 2 + V 4 ). 

Since equation 2 and 3 don't transform into each other as a result of the trans- 
position V 2 -H- V 3 [obviously, their coincidence with any other equations (jl]) 
is excluded], this transformation cannot belong to the group Go- 

In this paper, we don't discuss any methods for revealing the symmetry 
group G of the given dynamical model OH) (one such method was devel- 



oped in [16|). Let us only note that all symmetry elements of the system 
equilibrium state must enter the group Go. 

Let us consider a dynamical regime in our electrical chain @, described 
by the vector V(t) = {Vx(t), V 2 (t), ..V n (t)}. Acting on this vector by all 
elements g G Go, we choose only those which don't change it, i.e. elements 
satisfying condition = V. Full set of such elements represents a subgroup 
G of the group Go'. G C Go- This subgroup is the symmetry group of the 
considered dynamical regime V(t). 



According to the theorem proved in [18| , the variational system obtained 
by linearization of nonlinear equations (jl]) in the vicinity of the dynamical 
regime V(t) is invariant with respect to the above subgroup G C Go. This 
means that some elements g G Go which conserve the original system (jl]) 
don't conserve the corresponding variational system. 

The action of a symmetry element g G G on the vector V can be replaced 
by the action of a certain N x N square matrix. The set of such matrices 
associated with all elements g G G, acting in iV-dimensional space of state 
vectors V, represents a matrix representation T of the group G. 

The condition that G is the symmetry group of the variational system 
8 = J(t)S (this is the form of the variational system for the FPU-/3 model) 
leads to commutation of the matrix J(t) with all matrices of the representa- 
tion T [18]. This fact allows one to apply the well-known Wigner theorem 
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describing the structure of the matrix commuting with all matrices of a re- 
ducible representation of a given finite group. 

The conclusion of the Wigner theorem can be explain as follows. Let 
us decompose the reducible representation Y of the group G into irreducible 
representations (irreps) Tj of the same group: 



£%r,- (2i) 



(rrij are nonnegative integers). Such decomposition can be fulfilled by the 
appropriate transformation of the basis of iV-dimensional vector space corre- 
sponding to T (carrier space). Then matrix J(t), commuting with all matrices 
of T, acquires a block-diagonal form as a result of the above-mentioned trans- 
formation of the original basis. Every block Dj corresponds to a certain irrep 
Tj and possesses dimension equal to rij-rrij, where rij is the dimension of Fj, 
while rrij [see Eq. (I2~lj) ] determines how many times Tj enters the reducible 
representation F. Moreover, the block Dj represents a matrix with very spe- 
cific structure: it is a direct product of a certain rrij- dimensional matrix and 
n^-dimensional identical matrix. 



As was shown in |18|], this structure of the matrix J = J(t), obtained after 
basis transformation producing the form (121]) of the representation T, leads 
to splitting (decomposition) of the original variational system 8 = J(t)S into 
a number of independent subsystems whose dimensions can be considerably 
smaller than that of the original variational system. Namely, rij identical 
subsystems of dimension rrij correspond to every block Dj and, therefore, 
to the irrep Tj. Note, that induction numbers rrij entering Eq. (I2ip can be 
found with the aid of the following formula 

m ^ = AnJ2^9)XrX9) (22) 

11 11 geG 



from the theory of finite groups [31] . Here, ||G|| is the order of the group G 
(the number of its elements), while Xr{g) and Xr 4 (fi0 are traces of matrices 
corresponding to the element g in the representations F and Tj, respectively. 
The bar over Xr^d) denotes complex conjugation. 

One can obtain the explicit form of variational system 6 = J(t)5 decom- 
posed into independent subsystems by passing from the original basis in the 
space of all infinitesimal vectors 6 to the new one which represents the set 
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R of all basis vectors of the irreps Tj entering the decomposition (I2ip . This 
transformation of the basis leads to transformation of the old infinitesimal 
variables 5i{t) to the new one yi(t) as follows: 

V = S6, (23) 

where the columns of the orthogonal matrix S are vectors from R. 

Let us illustrate the above described method with a simple example. 

For the case N = 4, the variational system of the electrical chain (14"1) 
reads: 





(1 - V?)5\ - 


-2{V?)8\ 


= 2{V 1 V 1 


- l)6i + 5 2 + S 4 , 


< 


(1 - V 2 2 )8 2 - 


-2(V 2 2 )5 2 


= 2(V 2 V 2 


- l)5 2 + 5 3 + (5i, 




(1 - vi% - 


- 2(^ 2 )4 


= 2(VsV 3 


- 1)5 3 + 5 4 + 5 2 , 




k (i - vl% - 


- 2(Vi)5\ 


= 2(V A V A 


- 1)64 + 61 + 5 3 . 


EereV i = Vi(t),8 i = 6i(t),i 


= 1..4. 






Let a, = 1 




-2(^ 2 ),7, 


= 2{V i V i 


— 1). Then Eqs. 



written in the matrix form: 

where matrices a(t), /3(t) and j(t) read 



&(t) 



(a x 

a 2 



\0 






03 







a 4 J 



(Pi 








o\ 



















A» 





\o 











(24) 



can be 



(25) 



(26) 



7i (t) 



l(t)-- 
[li 
72 
73 
\0 



+72, 










A) 


1 





i\ 




1 





1 





,72 = 





1 





1 


/ 


V 





1 


0/ 



One possible dynamical regime in the model (jlj) is 

V(t) = {^i(t), u 2 (t),u 2 {t), U^t)}. 



(27) 
(28) 



(29) 
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This fact can be checked by straightforward substitution of this vector into 
Eqs. fll]). The four equations with respect to variables Vi(t),i = 1..4 reduce 
to two equations for new variables Ui(t) and U 2 {t). 

Note that the vector fl29|) determines a quasiperiodic dynamical regime, 
and we want to emphasize once more that our method can be applied not 
only to periodic, but also to quasiperiodic regimes. 

It is easy to find the symmetry group G of the dynamical regime (|21?|) . 
This is a second-order group 

G = [e, i] (30) 
whose nontrivial element % acts on the arbitrary state vector V as follows: 

i{Vi(t), V 2 (t), V 3 (t), V 4 (t)} = {V A (t), V 3 (t), V 2 (t), V l (t)}. (31) 

It transposes voltages as V 2 -f* V3, V± V4 and turns out to be a certain 
symmetry transformation of the model (jlj). Obviously, this transformation 
is simultaneously the symmetry element of the regime (j29|) . 

As any second-order group, the group (1301) possesses two irreducible rep- 
resentations, r x and T 2 , both one-dimensional, which we present in Table [TJ 



Table 1: Irreducible representations of the group G = [e,i] 



Irreps 


e 


i 




1 


1 


r 2 


1 


-1 



Now we must construct basis vectors, cp and if), of these irreps, which 
represent certain four-dimensional vectors. Let us remember the definition 
of the representation T (reducible or irreducible) of an arbitrary finite group 
G: 

g<(> = M{g)<t>, VgeG (32) 

Here 4> = {<j>ii 4>2-> ••■> ^n} * s ^ ne "supervector", representing the list of basis 
vectors of a given n-dimensional representation, M(g) is the matrix asso- 
ciated in this representation with element g G G, while M(g) denotes the 
transpose of the matrix M(g). Eq. f|32|) represents a system of linear alge- 
braic equations for unknown vectors 4 > i ,i = l..n. 

Traditionally, basis vectors of an irreducible representation are found by 
means of projection operators (see, for example [3JJ). However, it is more 
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convenient for our purpose to obtain these vectors using the straightforward 
method [32J based on the definition (I32|) of matrix representation (some de- 
tails of this method can be found in |18j). Obviously, it is sufficient to use 
equations fl32|) only for the generators of the group G. 

Let us find the basis vectors and i/> of the irreps Ti and T2, whose 
one- dimensional matrices are given in Table 1. 

From Eqs. ( 13"2|) we obtain for Fx and T 2 the following equations 



r x : iip = (l)<p, (33) 

T 2 : = (-l)t/>. (34) 

We search cp and in the form cp = (Vx, V 2 , V 3 , Vx) and (Ux, U 2 , U 3 , C/ 4 ), 
where Vi, U%{i = l.A) are arbitrary scalar values. From Eq. (|33|) . we find 

itp = z(Vx, V 2 , V 3 , V,) = (V 4 , V 3 , V 2 , Vx) = (Vx, V 2 , V 3 , V 4 ). 

Therefore, Vx — V4, V 2 — V 3 and 

cp = (Vx,V 2 ,V 2 ,Vx). (35) 

From Eq. flMj) . we obtain 

itP = i(Ux, U 2 , U 3 , Ux) = (U 4 , U 3 , U 2 , Ux) = (-Ux, -U 2 , -U 3 , -Ux). 

Therefore, Ux = —Ux, U 3 = —U 2 and 

*l> = (Ux,U 2 ,-U 2 ,-Ux). (36) 

In our example, both vectors if) and ip represent two-dimensional subspaces 
of the space of all state vectors, and they can be written in the form: 



<p = V 1 



i> = Vx 






VJ 




v-v 



v 2 



+ v 2 



M 
1 
1 

1 

-1 



(37) 



(3? 



13 



Each of these vectors determines two basis vectors of the four- dimensional 
space of the nonlinear model (jl]) or linear model OH]) for N = 4. Normalizing 
these basis vectors, we can write them as a columns of the matrix S entering 
Eq. ([23]): 


1 
1 




S = —= 



1 
3 



ft 







V 



i 






(39) 



0\ 
1 

-1 
■1 0/ 

The columns of the matrix S represent the orthogonal and normalized basis 
in the space of all variables of the dynamical models (j3J) and ( |24"|) for the case 
N = 4. 

Using the matrix (1391) . we can split the variational system in the form (J2! 
or ( 1251) into two independent subsystems. 

Indeed, for dynamical regime ( )29|) . matrices a, /3, 7 acquire the forms 
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(40) 
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7(*) 



(45) 



/71 + I 1 \ 

1 72 + I 

7! - 1 1 

1 72-I/ 

Each of the above matrices can be written as direct sums of certain two- 
dimensional blocks. Therefore, the variational system ([21]) is decomposed [in 
the variables i/j (J = 1..4) from Eq. into two independent subsystems 
of the form: 



cti(t) 






a 2 (t) 

ai(f) 
a 2 {t) 



y 



y + 














y 



y 



7,(i ! +1 



where y = (2/1,2/2)- 

Finally, let us note that as a result of the relations ( 1371) and (I38j) . one can 
write the decomposition formula ff2T]) for our case as follows: 



r = 2Ti + 2r 2 



(47) 



4. Variational system splitting for the NNMs ffT3l)- flTTj) 

Each nonlinear normal mode from Eqs. (jl3p - (ll7p possesses a certain sym- 
metry group G which contains a translational subgroup T n . Let us consider 
the mode (f> 2 . It can be excited in chains for which n = ^ is an integer 
number (N mod 4 = 0). The mode (f) 2 [see Eq. (JHJ)] is determined by ex- 
tended primitive cell (EPS) of the form {V(t), V(t), -V(t), -V(t)} which is 
repeated along the chain. This EPS being four times larger than the cell 
in equilibrium state describes a specific vibrational regime of voltages in the 
electrical chain (J4]) under the periodic boundary conditions. 

Obviously, translation of the voltage distribution (fT4"|) by 4a, where a is 
the cell size in equilibrium state, does not change the form of this distribution 
and, therefore, a 4 is one of its symmetry element. This element is a generator 
of the translational symmetry T n . 

The same translational group corresponds to the mode (f) 3 f fl~5|) . However, 
the full symmetry groups of the modes 4> 2 and 3 are different: they differ 
by the second generator, iu for (p 2 and ai for </> 3 . 

One can split variational system for (f) 2 and (f) 3 by the full symmetry 
groups of the corresponding dynamical regimes, [a 4 ,m] and [a 4 , cm], but for 



15 



our purpose it is more convenient to use a common group G = [a 2 u] of these 
modes which is defined by the generator a 2 u. The action of this symmetry el- 
ement on the voltage distribution {Vj(t)\j = 1..N} reduces to 2a-translation 
accompanying by change of sign of all variables Vj(t). Note that square of 
the symmetry element a 2 u is equal to the above discussed translational gen- 
erator a 4 . Using the group G = [a 2 u], we can split the variational systems 
for both NNMs (fi 2 and </> 3 with one and the same matrix S. 

Now according to the general prescription, we must construct the ba- 
sis vectors of the irreps of the group G. All such representations are one- 
dimensional, because this group is cyclic and, therefore, abelian. Every irrep 
Tj of G can be fully determined by the one- dimensional matrix jj associated 
with the generator a 2 u of this group. The number of irreps of the group G is 
equal to its order n — y. We can define these irreps, Tj, by 1 x 1 matrices: 

li = e_ ™~; 3 — 0, 1, 2, ...{n — 1) [jj are roots of unity order n\. 

As was described in the previous section, the basis vectors ?/> ■ of one- 
dimensional irreps 1^ can be obtained by solving the following equation 

a 2 ui^ j = -fjifij, (48) 

where ij)^ = {V\, V2, ..Vat} is an arbitrary voltage distribution for our elec- 
trical chain. For the considered case, we can easily find (see some details 
in 18]) that the solution to Eq. fj48|) represents a two-dimensional space L 



of all possible voltage distributions which can be written in the form 

| 7 7 4 a, 7 - 4 &,-7- 5 a,-77 5 &l-}, 1 ' 

where a and b are arbitrary constants. Setting a = 1, b = and a = 0, b = 1, 
we select two basis vectors of the above subspace. 

Finally, looking over all values j = 0, 1, 2, (y — 1), we find a certain 
basis of the full space L: 



<tl>f = {|1, 0, - 7 -\ 0| 7 7 2 , 0, - 7 " 3 , 0|7" 4 , 0, - 7 " 5 , 0|...}, 
*/>f = {|0, 1, 0, -t^IO, lj\ 0, -77 s , 0|0, 7 " 4 , 0, -77 5 |-} 



(50) 



After normalization, we can use vectors (|50|) as columns of the matrix 
S which splits the variational systems for 2 an d 03 into independent sub- 
systems. Note that taking into account any additional symmetry elements 
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beyond those of the group G = [d 2 u] can generally reduce dimensions of the 
above subsystems that were obtained with the aid of the group G. 

Proceeding in such a manner, we reveal that the variational systems for 
studying stability of NNMs f|T3|) - f|T7|) can be split into independent subsys- 
tems whose dimensions are equal to 1, 2, 2, 3, 3, respectively 

Below, we present these subsystems in matrix form 

A8 + B8 + D6 = 0, (51) 
where matrices B, C, D are constructed from the following elements: 

c(t) = \-V\ 

f(t) = 2(V 2 + VV-l), (52) 

g(t) = -4VV. 

Here V = V(t) is the solution of the governing equation (TSUI) for the corre- 
sponding nonlinear normal mode. 

1) = v-(t){i, -1|1, -1|1, 

For this case all subsystems f[5"Tj) are one-dimensional with the matrices 
A = c(t), B = g(t), D = f(t) - 2 cos(A;), (53) 

where k — 

For NNMs 2 and 3 , subsystems (|5ip are two-dimensional: 

2) 2 = V(t){l, 1, -1, -1|1, 1, -1, 

HWO-Hvr)- (54) 

3) 3 = y(t){i,o,-i,o|i,o,-i,o|...}. 

H !)■*-(! S)-HV^)- w 

In Eqs. (JMD-dSSD, 7 = e^ j . 

4) For both NNMs 4 = V(t){l, 0, -1|1, 0, -1|...} and 

<p 5 = V(t){l, 1, 0, — 1, — 1, 1 1 , 1, 0, — 1, — 1, 0} three-dimensional subsys- 
tems flBTl) turn out to be identical. These subsystems are determined by the 
matrices 

(c 0\ (g 0\ // 1 7 \ 

i= 1 ,5= \,D= 1 -2 1, (56) 
\0 cj \0 \7 1 // 
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where 7 = e^~K 

Thus, stability of all NNMs in the electrical chain (jlj) can be investigated 
by analysing a set of one-, two- and three-dimensional systems of differential 
equations with time-periodic coefficients determined by the function V(t). 

5. Stability analysis of 7r-mode 

As we have just seen, the most simple splitting takes place for 7r-mode. 
Indeed, in this case the original variational system for the chain with any 
everU number of cells (N) can be split into one- dimensional subsystems which 
are determined by Eqs. (l5Tj) -( l53"]) . In the explicit form these equations read 

(1 - V 2 )5 - A(VV)5 - 2(V 2 + 2VV)5 = -4 cos 2 (|)5, (57) 

wherefc = Mj = l,2,...,f. 

Eq. (15T|) is a linear second-order differential equation with first derivative. 
One can reduce it to the so-called normal form which doesn't contain the first 



derivative by the appropriate changing of variables S and t [33]. In our case, 
this procedure is simplifyed essentially, because the left part of Eq. ( 1571) may 
be written as the second derivative of the expression [(1 — V 2 )<5]. Therefore, 
introducing the new infinitesimal variable 

y(t) = [l-V 2 (t)]5(t), (58) 

we immediately obtain a very simple equation 

y + 4c(t)cos 2 (|)i/ = 0, (59) 

where c(t) = [1 - V 2 (t)]-\ 

Let us comment on the notation used in this equation. The time-periodic 
function 

c(t) = l/[l-V 2 (t)] (60) 
is determined by governing equation (1201) with fi = 4 for initial conditions 

V(0) = A, V(t) = 0. (61) 



2 Note, that condition iV mod 2 = is necessary for existence of 7r-modc. 
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Thus, A is the amplitude of voltage oscillations on capacitors. The variable 
k = jfj, j = 1..N is introduced instead of the mode number j (for the limit 
case iV — > oo it represents the one-dimensional wave vector). 

As was shown in Ref . 1(| , the governing equation (120]) describes a periodic 



dynamical regime for A < 1. The value A = 1 represents singular point 
beyond which the motion becomes infinite. 

In FigfJl we present plots of function V{t) for three values of its amplitude: 
A = 0.1, A = 0.8, A = 0.99. As one can see, practically harmonic oscillations 
take place for small amplitude (A = 0.1), while nonlinearity of vibrations 
is observed clearly for A = 0.8. For A — > 1 oscillations acquire specific 
triangular form. 

The physical cause of the stability loss of 7r-mode, as well as other NNMs, 
is parametric resonance. Let us consider this cause for the case of small 
vibrational amplitudes. For A — > one can neglect all nonlinear terms in 
Eq. fTSOj) and this equation then acquires the form 

V + W = 0. (62) 

This is equation of the harmonic oscillator with eigenfrequency u = 2 and, 
therefore, we can write its general solution in the form: 

V(t) = Asin(2t + ^ ), (63) 

where po is an initial phase. 

As we have already stated, the variational system for the Floquet analysis 
of 7r-mode can be decomposed into the set of independent equations (1591 
with different values of the parameter k. As a consequence of the condition 
of small amplitudes, \V(t)\ <C 1 and one can decompose the function c(t) = 
[1 — V 2 ^)]^ 1 into Taylor series. Taking into account only first nonlinear term, 
we obtain 

5 + A 1 [l + V 2 (t)]5 = 0, (64) 

where 7 = cos 2 (|). Introducing new time variable r = 2t + ip , allow us to 
transform this equation to the standard form of the Mathieu equation jjil ]: 

y" + {a-2q cos(2r)}y = 0. (65) 

Here 

Vi r ) = $ C 2 ° ), a = (1 + -yh> q = (66) 



19 



while differentiation with respect to the new time argument is denoted by 
prime. 

It follows from Eq. f )66|) that there is a linear relation between parameters 
a and q of the Mathieu equation: 

a = 7 + 2q. (67) 

On the other hand, it is well known that there exists an infinite number 
of regions of unstable motion in the (a — q) plane for the Mathieu equa- 
tion (IBISj) (see Ince-Strutt diagram in [34J]), and these regions appear because 
of parametric resonance. 

For studying stability loss of 7r-mode, it is sufficient to consider only the 
first region of instability from the Ince-Strutt diagram. We depict it by 
grey color in FigfSl We also depict in this figure a number of straight lines 
a = 7 + 2q [see Eq. (1671) ] for different values of the parameter 7, which is 
connected with the mode number j by the relation 

7 = cos 2 (^). (68) 

For A = 0, the function c(t) from Eq. (|60|) turns out to be equal to unity 
and equation ( l59l) represents a set of harmonic oscillators with eigenfrequen- 
cies 

Uj = cos(|) = cos(^), j = 1..N. (69) 

On the other hand, for A 7^ 0, Eq. ( 1571) transforms to a set of the Mathieu 
equations with parameters a and q depending on the mode number j. 

The function c(t) is time-periodic and for certain mode numbers j the 
coefficients a(j) and q(j) may hit the unstable region depicted in Figj3j In 
such a case condition of the parametric resonance is satisfyed and, as a result, 
the corresponding 5j(t) begin to grow exponentially. This means the loss of 
7r-mode stability, because this mode can be stable only when all Sj(t) continue 
to be small if their initial values Sj(0)(j = 1..N) are small. 

Let us consider the interaction of 7r-mode with a " sleeping" linear normal 
mode 5 jit), i.e. a mode equal to zero at the initial instant t — 0. If we 
gradually increase the 7r-mode amplitude A from the value A = 0, the point 
(a, q) moves upwards along one of the straight lines depicted in Fig IS] from its 
initial position (7, 0), because a = r )+2q and q = 7^ [see Eqs. (IB"?]) and (1661) ]. 
The 7r-mode is stable up to intersection of the above straight line with lower 
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boundary of the grey region in FigJ3J This intersection corresponds to the 
maximal amplitude A 3 C of 7r-mode for which excitation of the sleeping mode 
5j(t) still does not occur. 

On the other hand, the value A 3 C can be obtained as a solution of the 
equations 

a = 1 -^¥-' (70) 
a = 7 + 2q, 

where the first equation defines the lower boundary of the gray region in 



Figj3] (see, for example, |34j). As a result of solving Eqs. ( I7D1 . we fin< 



J3 



« = i^ (7i) 

Taking into account the relation q = 7^ and Eq. we obtain finally 

Ai = ^4). (72) 

Obviously, the loss of 7r-mode stability is determined by that sleeping 
mode 5j(t) for which minimal value of A J C (j = 1..N) corresponds. Then we 
conclude from Eq. ( I72|) that this minimal value A c corresponds to j = 1, 
i.e. the mode 5\(t) is excited firstly when 7r-mode amplitude A begins to 
increases from zero. Hereafter, we call A c critical 7r-mode amplitude because 
this mode loses its stability beyond A c . It follows from Eq. (172p that for 
7r-mode A c — > when N — > 00. 

In FigJU we depict for 7r-mode the plot of the function A C (N) which 
shows dependence of the critical amplitude A c for different length (N) of the 
considered electrical chain. Solid line in FigJU corresponds to the function 
A C (N) obtained from Eq. ( 172|) for j = 1, while dotted line is depicted with 
the aid of the 7r-mode Floquet stability analysis. We calculate the Floquet 
exponents integrating Eq. ([57]) for a given amplitude V(0) = A and for all 
possible values of 7 detecting instability of 7r-mode if at least one of these 
exponents exceeds unity more than e = 10~ 5 . 

Both plots in FigJU are in very good agreement for large values N. Some 
deviations taking place for small iV can be explained by approximation that 
was done for simplification of the governing equation (|20p . 



5 For \q\ -C 1 one can neglect the term %- in Eq. (|70|) . 
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For j = 1 and iV ^> 1 , we find from Eq. (I72p 



Therefore, decreasing of A C (N) for iV — > oo is determined by power the law 

A C {N) = CN~P, (74) 

with C « 3.629, /3 = -1. 

Obviously, the above 7r-mode stability analysis based on the transition 
from Eq. (l5Tj) to the Mathieu equation turns out to be exact in the thermo- 
dynamic limit N — > oo. 



6. Stability analysis of the nonlinear normal mode <p 2 

According to Eqs. fl5Tj) and fl54|) . stability of NNM <p 2 is determined by 
stability of zero solution of the two-dimensional systems with different values 
j = e * J . 

However, it is easy to note that matrices A,B,D from Eq. f lo^|) commute 
with each other and, therefore, they can be diagonalized simultaneously with 
the aid of an appropriate unitary transformation. As a result of this trans- 
formation, we split each of these two-dimensional variational systems into 
two independent scalar differential equations. 

Indeed, Eqs. ( ]5Tj) and ( 154")) can be written in the form: 

[1-V 2 )5 1 -4VV6 1 -(2V 2 + VV)5 1 = -2^ + (1 + 7 )5 2 , 
1 - V 2 % - WV5 2 - {2V 2 + VV)5 2 = -28 2 + (1 + j)S v 

Introducing new infinitesimal variables yi — (1 — K 2 )5j instead of the old 
variables S{ (i = 1,2) [see Eq. (!23|) ]. we obtain 

1-^=^1 + 7 -2 J*' (76) 
i.e. 

My = Ry. (77) 
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With the aid of a certain unitary transformation T, one can diagonalize 
matrix R without changing matrix M which is proportional to the identical 
matrix. As a result of such transformation, we obtain: 



1 -V 2 \ .. (Ai 



1 - V 2 V A 



\2 



e, (78) 



where e = Ty, while Ai = — 2 + 2sin(|), A 2 = — 2 — 2sin(~) are eigenvalues 
of the matrix R. 

Thus, stability analysis of the mode (f) 2 is reduced to studying stability 
of zero solutions of the independent equations: 

ei + 2c(*)[l-sin(-)]ei = 0, 

k (79) 
£ 2 + 2c(t)[l + sin(-)]e 2 = 0, 

where c(t) = (1 — V 2 )^ 1 . It is obvious, that these equations can be investi- 
gated in the thermodynamic limit iV — > oo in the same way as the variational 
equation (}59l) for 7r-mode (note, that k = j^j for 7r-mode, while in the case 
of </> 2 -mode, k = -jfj). The corresponding result read: 

The scaling of A C (N) in the limit N — > oo for the mode (f> 2 , as well as that for 
the mode is determined by Eq. ( 174)) . However, in this case, the exponent 
j3 is equal |. The final stability result for NNM (f> 2 is depicted in Figj5j 

We have reached the complete splitting of the variational system for NNM 
2 in two steps. Firstly, we have used the general group-theoretical approach 
based on the symmetry group G = [a 2 u] and obtained two-dimensional sys- 
tem (Ff6l) . Secondly, we have split these system into individual scalar equa- 
tions (I79p . Let us note that one can obtain this final result in one step only, 
if splitting is produced on the base of the full symmetry group G 2 = [a 4 ,iu] 
of the mode (f> 2 . 

7. Stability analysis of NNM 3 

Stability analysis of nonlinear normal modes 3 is reduced to studying 
two-dimensional system, determined by Eqs. flBT]) . (l55p. which cannot be split 
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into scalar equations. Indeed, matrices A and B don't commute with the 
matrix D : AD ^ DA, BD ^ DB. As a consequent, it is impossible to 
diagonalize these three matrices, A, B, D, simultaneously, i.e. with the aid of 
one and the same unitary transformation, as it was possible in the case of the 
mode 02- m turn, this means that two-dimensional variational subsystem 
for NNM cf> 3 cannot be decomposed into independent scalar equations and 
we are forced to study zero solution stability for the following 2D system: 

c{t)8\ + g(t)5i - - (1 + 7 )5 2 = 0, 

<5 2 -25 2 - (1+7)61 = 0, 

where 7 = e tk {k = -jf)- The time- dependent functions c(t),g(t) and f(t) are 
defined in Eq. fl52|) . 

We can consider infinitesimal variables 81 and 82 as linear normal modes 
(phonon modes) corresponding to certain wave numbers. Let us remind, that, 
rigorously speaking, every phonon mode represents time-dependent variable 
5k{t) multiplied by a certain iV-dimensional vector determining a distribution 
of voltages on all capacitors of the electrical chain. We call Sk{t) "phonon 
mode" only for brevity. It follows from the matrix S, which leads to splitting 
of the original variational systems, that the above normal modes, 81 and 82, 
are associated with different wave numbers, k and tt — k: 

81 = 8 k (t), 8 2 = Sir-kit). 

It is interesting to note that the wave numbers k and tt — k situated 
symmetrically with respect to the center of primitive cell [0; tt] of the recip- 
rocal lattice which corresponds to the original lattice of our electrical chain 
in vibrational state describing by the mode (f>$2. 

As in the case of 1 - and 2 " m °des, NNM <f> 3 turns out to be stable for 
sufficiently small amplitudes A = V(0). Increasing of the amplitude leads 
to parametric resonance appearance which results in exciting two sleeping 
phonon modes Si = 5k{t) and 82 = 8 n -k(t) simultaneously. 

The stability loss of NNM 3 is determined by those sleeping phonon 
modes (the modes with certain numbers j) which are firstly excited as a result 
of increasing the amplitude of this nonlinear normal mode. Calculating the 



4 The primitive cell [0; — ] in the reciprocal lattice corresponds to the one-dimensional 
lattice with period a. 
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Floquet exponents, for electrical chains for different N, we can detect the 
corresponding phonon modes. In turn, this give us a possibility to obtain the 
function A C (N) similar to that depicted in FigJUfor the mode cf> 1 . 

8. Stability analysis of modes <fi 4 and <p 5 

The stability analysis of nonlinear normal modes 4 and d> 5 is reduced 
to studying zero solution stability of the one and the same three-dimensional 
system determined by Eqs. f joTj) and ( 151)]) : 



( c{t)5 3 + g{t)5 3 - 7ft - 5 2 - S 3 = 0. 
Here the functions c(t), g(t), f(t) are given in Eq. fl52l) . while 7 = e lk with 



Results on stability of NNMs </> 4 , <f> 5 are presented in FigJBJ It is inter- 
esting to note that the stability plots for the modes (f) 4 and 5 turn out to 
be identical, despite the different functions V(t), determining the coefficients 
c(t),g(t), fit), correspond to these modes. Indeed, the functions V{t) for the 
modes 4 and 4> 5 are obtained from the governing equation (|2"01) with ji = 3 
and // = 1, respectively. 

9. Asymptotic behavior of A c (iV) for N — >• 00 

The critical amplitudes A C (N) for all symmetry-determined nonlinear 
normal modes (fTBl - ffTTj) tend to zero when iV — > 00. For the mode <p l and <f> 2 
the corresponding scaling law of decreasing A C (N) was presented in Eq. (1711) 
with (3 equal to 1 and |, respectively. Naturally, it is interesting to find simi- 
lar asymptotic formulas for other NNMs. Numerical experiments allow us to 
obtain Figj6] where we present functions A C (N) in logarithmic scale for large 
N (N > 100). It is obvious from the plots depicted in this figure that A C (N) 
for all NNMs are power functions since dependence of In A C (N) on In iV rep- 
resents straight lines. The coefficient /3 in the formula A C (N) = CN 13 is 
determined by inclination of the corresponding straight line to the horizontal 
axis. 




(82) 



fc = ^(j = l,2,...,f) 
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In this way, we have obtained the following values of the constants (3 and 
C entering the law A C (N) = CN 13 : 



1 [a 2 ,m] — vr-mode [Eq. (1T3"]) ]: = 0.99, d = 3.62; 
2 [a 4 , z«] [Eq. (ED]: /3 2 = 0.49, C 2 = 2.89; 

03 [a 4 , at] [Eq. ([15])]: /3 3 = 0.99, C 3 = 8.57; 1 J 

4 [o 3 ,zu], 5 [a 3 w,aw] [Eqs. (JTSJ and (pi])] : /3 4 = 0.99, C 4 = 7.54. 

Relying on the above data, we can assert that for all NNMs, except for 
the mode 2 , the constant /3, determining the rate of decreasing of the criti- 
cal amplitude A c with increasing N, is equal to unity ((3 = 1), while for 2 
it is equal to |. Note that we have obtained the values of the parameter 
f3 for the modes </> 1 and 2 analytically, as well as numerically, while for 
NNMs 3 , 4 , 5 this parameter was found only numerically. We conjecture 
that it can be obtained also analytically by the specific asymptotic method 



developed in |35j. However, this was not done in the present paper, because 
numerical results ( 183]) seem to be sufficient for hypothesis that the parameter 
(3 for NNMs 3 , 4 and 5 is equal to unity. Indeed, the numerical data ( 183]) 
show that the deviation of the parameter (3 from 1 for modes 1; 3 , 4 , 5 
and that from | for the mode 2 are equal to 0.01. We associate these devi- 
ations with some computational errors, because numerical data demonstrate 
the same deviation from the exact values j3\ = 1, (3 2 = ~ which were found 
in Sees. 5, 6. Naturally, we can assume that exact values of 03,134 and [3$ are 
equal to unity. 



10. Conclusion 

In the present paper, we study periodic vibrations in the chain of nonlin- 
ear capacitors coupled by linear inductors, assuming that volt age- dependence 
of capacitors is described by the function C(V) — C (l — bV 2 ) as was done in 
Ref. [lo[. These dynamical regimes represent symmetry-determined nonlin- 
ear normal modes and we prove that only five types of such modes can exist 
in the considered electrical chain. The stability analysis of these nonlinear 
modes with the aid of the straightforward Floquet method turns out to be 
very difficult in the case of large number (N) of the electrical chain cells. 

We apply the general group-theoretical method, developed in for 
a radical simplification of the stability studying. This method allows us 
to reduce large-dimensional variational systems, appearing in the stability 



26 



analysis of nonlinear normal modes, to independent subsystems of small di- 
mensions. Even in the case N — > oo for all such modes it is sufficient to 
study zero solution stability of one-, two-, and three-dimensional variational 
systems. Using this approach, we were able to find the voltage amplitude 
A C (N) of a given NNM which corresponds to onset of unstable vibrations 
in the electrical chain for any fixed N. Moreover, we justify the scaling law 
A C (N) = CN 13 with parameter /3 equal to | for the mode 2 [ see Eq. (|14p ] 
and equal to unity for all other NNMs. 

It is essential that our group-theoretical method can be applied not only 
to stability studying of periodic dynamical regimes, but also to other regimes. 
In particular, this is true for bushes of nonlinear normal modes which repre- 



sent quasiperiodic vibrations [25|, |27J, |28[ . In such a case, one cannot apply 



the standard Floquet method. However the variational systems correspond- 
ing to these dynamical objects can be also decomposed into subsystems of 
sufficiently small dimensions and this simplifies the stability analysis. For 
the case of electrical chain (|4j), this problem will be considered in a future 
paper. 
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Figure 1: Chain of nonlinear capacitors coupled by linear inductors 
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Figure 2: Function V(t) for three values of 7r-mode amplitude: a) A = 0.1; b) A = 0.8; c) 
A = 0.999 




Figure 3: First zone of instability for the Mathieu equation and a number of straight lines, 
determined by Eq. (1571) 
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Figure 4: Dependence of the critical amplitude A c on the number TV of the chain cells for 
7r-mode. Solid line corresponds to the function A c (t) obtained from Eq. (|72|) for j = 1, 
while dot-dashed line was obtained as a result of numerical calculation 
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Figure 5: Stability of nonlinear normal 
modes (|13p - (fl7)) for electrical chains 
with different N. Solid line corresponds 
to the mode 2 CHJ , dashed line - to 



the mode </> 3 ([T5]) . dot-dashed line - to 
the 7r-mode <p ± (fT3"]) , dotted line - to the 
modes 4 (|T6|) and 5 (fT7|) 
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Figure 6: Plots of A C (N) for nonlin- 
ear normal modes (|13[) - (|17[) in logarith- 
mic scale. Solid line corresponds to the 
mode 4> 2 $M , dashed line - to the mode 
4>3 (|15[) , dot-dashed line - to the 7r-mode 
4>i C31), dotted line - to the modes <p 4 
dTCJ) and 5 (HID 
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